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The capacitance of mesoscopic samples depends on their geometry and physical properties, de- 
scribed in terms of characteristic time scales. The resulting ac admittance shows sample to sample 
fluctuations. Their distribution is studied here -through a random-matrix model- for a chaotic cav- 
ity capacitively coupled to a backgate: it is obtained from the distribution of scattering time delays 
for the cavity, which is found analytically for the orthogonal, unitary and symplectic universality 
classes, one mode in the lead connecting the cavity to the reservoir and no direct scattering. The 
results agree with numerical simulations. 



PACS numbers: 72.20.My, 05.45.+b, 72.15.Gd 

The elementary notion of capacitance of a system of 
conductors, as a quantity determined solely by the ge- 
ometry, has to be revised if the electric field is not com- 
pletely screened at the surface of the conductors. In fact, 
the penetration distance of the field is of the order of the 
Thomas-Fermi screening length, which may be apprecia- 
ble for a mesoscopic conductor: the standard description 
of a capacitor in terms of the geometric capacitance C e 
(that relates the charge Q on the plate to the voltage U 
across the capacitor), gives way, in the mesoscopic do- 
main, to a more complex entity C^, the electrochemical 
capacitance (that relates Q to the electrochemical poten- 
tial of the reservoirs), which depends on the properties of 
the conductors JlJ. This fact, in turn, has important con- 
sequences for the ac current induced in the system when 
the electrochemical potentials are subject to a nonzero- 
frequency time variation . 

The electrochemical nature of the capacitance has been 
relevant to a number of experiments Q and has been dis- 
cussed theoretically by several authors Remark- 
ably, it has been found that the resulting ac admittance 
can be described in terms of characteristic time scales re- 
lated to energy derivatives of scattering matrix elements. 

It is well known that, as a result of quantum interfer- 
ence, the dc conductance of mesoscopic structures shows 
strong fluctuations as a function of the Fermi energy or 
the magnetic field, as well as from sample to sample. A 
statistical analysis of this phenomenon has been done, 
for diffusive transport in disordered structures, using mi- 
croscopic perturbative and macroscopic random-matrix 
theories JbJ , and for ballistic microstructures -cavities in 
which impurity scattering can be neglected so that only 
scattering from the boundaries is important- whose clas- 
sical dynamics is chaotic, using semiclassical, field theo- 
retic and random-matrix approaches . 

An extension of the above random-matrix studies to 
include the ac admittance of mesoscopic structures is the 



subject of the present investigation. 

In this letter we shall confine our discussion to the 
geometry shown in Fig. 1. In this system there is, of 
course, no dc transport, but there may be an ac current, 
determined by the admittance jjJU 



i + 



(1) 



written in the Thomas-Fermi approximation and to low- 
est order in the frequency uj. Here, g(ui) , g 1 (uj) de- 
note the admittance for the noninteracting and interact- 
ing system, respectively, the former being given, for zero 
temperature, by 
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Here, S(E) is the N x N scattering matrix for the system 
formed by the cavity and the lead, N being the number 
of propagating modes, or open channels, in the lead; A 
is the mean level spacing for the cavity (the inverse of 
the level density) . Following || , we have introduced the 
dimensionless time delay 
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where exp(i6 
as 



2ttA dE ' 

= det S . We then write g 1 (uj) of Eq. (| 



/(w) = -iuC e a H , (4) 

where the dimensionless capacitance a is given by 



a = CJC e 
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Notice that, for a macroscopic cavity, t| < 1, so that 
a ss 1 and ~ —iujC e . 

The one-energy statistical distribution of the S ma- 
trix for ballistic cavities larger than the Fermi wave- 
length has been modelled successfully through an "equal- 
a-priori probability" ansatz (known as a "circular ensem- 
ble") 0,0], when the classical dynamics is chaotic and 
direct processes through the microscructure can be ne- 
glected, so that, as a result, the averaged S vanishes, 
S = 0. It is clear, though, that the time delay r of Eq. 
(||) is a two-energy function and thus requires more in- 
formation for its statistical study. The distribution of 
t, w(t), has been studied for a one-dimensional disor- 
dered system within the invariant imbedding formalism 
in ||. In another approach, an underlying Hamiltonian 
described by a Gaussian ensemble was assumed and the 
problem analyzed using supersymmetry techniques: the 
two-point correlation function for the S matrix elements 
was derived in fiofl ; phaseshift times for unitary sym- 
metry, iV and S arbitrary were studied in Q. Ref. @ 
finds an approximation to w(t). We concentrate, in what 
follows, on w(t) for arbitrary symmetry (orthogonal, uni- 
tary and symplcctic, identified as (3 = 1, 2 and 4, respec- 
tively), N = 1 and S = 0: we show that this case can 
be treated using an old conjecture by Wigner Jl^,0 : we 
believe that the simplicity of the argument is appealing 
and gives an interesting perspective to the problem and 
a unified point of view for arbitrary (3. We also remark 
that, for ballistic cavities, the case of just one open chan- 
nel, N = 1, is very relevant from an experimental point 
of view, since cases of small N have been realized in the 
laboratory Il5|. We find below 
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where < r < oo. For (3 = 2, this result agrees with 
that of Ref. [pT| . The main result of the present paper, 
i.e. the (3 dependent distribution of the dimensionless 
capacitance a [a is related to the ac admittance via Eq. 
(0)1, then follows as 
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for 0< a < 1. A plot of pi^(a) for various values of r\ 
is presented in Fig. 2. For a macroscopic cavity, 77 — > 
and pa „{a) — > 8(1 — a). We now derive the distribution 
of time delays, Eq. (0). 
We write S for N = 1 as 
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For pure resonance scattering the K function can be 
given the sum-over-resonance form p^Jl4[] 
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where the "widths" T\ for a given symmetry class (3 can 
be written in terms of real amplitudes 7^ as 
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The quantity (E)/2 = h(E) was studied extensively 
by Wigner p3 14 1 ; it is called the "invariant derivative", 
because it remains invariant under the transformation 
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4> being a constant; since K = tan(6>/2), ( |12j ) takes 0/2 
to 6/2 + <p, and hence S to e^Se^. Both if and its 
transforms have the form K = tan J E h(E)dE, c being 
different for different transforms. Starting from one pole 
Ei of K , one can obtain the next one by determining 
the abscissa E2 so that the area under h(E) between 
Ei and E2 is w. Moreover, at a pole E\ we have T\ = 
l/h(E\). These relations are shown in Fig. 3. The levels 
and widths of the transforms of K can be obtained by a 
similar construction, starting at another abscissa. 
From (^,[n|) we find the energy average of S(E) as 



S(E) = S(E + il) = 
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where I -> 00 |6[ and t = nT/A. For S(E) = [circular 
ensemble, invariant under (113)] , we have t = 1. In this 
case (referred to in Ref. Jl^,[lJ] as that of a "normalized" 
R function) Wigner proposes the 

Conjecture: the statistical distributions of level spac- 
ing^ and residues are invariant under the transformation 

The above statemenent is a "conjecture" , not a "theo- 
rem" , and it is not clear, a priori, for what distributions, 
if any, it is fulfilled. Wigner, in his papers, proposes it for 
"most statistical distributions" . The conjecture, in rela- 
tion with the residue distribution, was verified numeri- 
cally for the case in which the energy levels entering Eq. 
( |Io| ) are constructed from a Gaussian Orthogonal, Uni- 
tary or Symplectic Ensemble, and the 7^ of Eq. ([ll]) as 
independent Gaussian variables: the residue distribution 
was found to remain invariant, within the statistical er- 
ror bars of the numerical simulation. On the other hand, 
the conjecture is seen, in our numerical studies, to be vi- 
olated for a spectrum of statistically independent energy 
levels following a Poisson distribution. 

Call Q (h) the probability density of the inverse widths 
h(E\) — h\ and P(h) the probability density of h 
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across the energy axis, irrespective of whether we are 
at resonance or not; P(h) is related to w(t) as w(t) = 
(tt/A)P(ttt/A). Assuming the above conjecture, Ref. 
H shows that 
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This relation can be understood by means of a very sim- 
ple argument. Consider, for one given K, the following 
level-average 
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for an arbitrary function /. From Fig. 3 we see that we 
cannot replace the sum in this equation by an integral. 
However, using the transformation ( |l2| ) we can construct 
"replicas" of K, all having the same distribution of h\; we 
do this n times, in such a way that the area between two 
successive levels is subdivided into n strips of area n/n 
each. Now we have a fine mesh, the sum over which can 
be approximated by an integral, using a density nh/ir, 
since the base of one of the above strips, at the place 
where h is the local value of the curve, is ir/nh. We then 
arrive at the above relation (|lj). 

If we use, in (14), the variable u = n/hA (and denote 



the distributions with a hat), we have 
P(u) = uQ(u) . 
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On the LHS, u can be thought of in terms of r of Eq. (||) 
as w = 1/r; at resonance, u takes the value u\ = ttT\/ A, 
which is the relevant variable on the RHS of Eq. (16). 
Thus, knowing the distribution of widths Q(u), Eq. ( |l6| ) 
allows finding P(u) Jrj| . 

For the three universality classes (3 = 1,2,4 and in- 
dependent Gaussian variables 7> , the distribution Q(u) 
is the chi-square distribution function with f3 degrees of 
freedom, 
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Eq. ( |16| ) then gives P(u), from which we find the dis- 
tribution of time delays u^(t) of Eq. ([?]). We notice 
the remarkable fact that, while wpir) certainly depends 
on the distribution of widths, other characteristics of the 
spectrum become lumped together in the invariance prop- 
erty contemplated in Wigner's conjecture. 

A numerical verification (using the simulation ex- 
plained above in relation with Wigner's conjecture) of 
w p( T ) of Eq. (0) is shown in Fig. 4 for the three symme- 
try classes (3 = 1,2,4: in all cases the agreement is seen 
to be very good. 

To summarize, we have found the statistical distribu- 
tion of capacitances pp^(a), Eq. (pft, a being defined in 



Eqs. (^J5j), for the system shown in Fig. 1, whose essen- 
tial element is a mesoscopic capacitor. The plate coupled 
to the backgate is a chaotic cavity; the experimentally 
relevant situation of one open channel (N = 1) is consid- 
ered and the possibility of direct reflection by the cavity 
is neglected. The essential ingredient that is needed is 
the statistical distribution wp{r), Eq. ([?]), of time de- 
lays r associated with the scattering from the cavity. It 
is shown that wp (r) can be obtained in a very simple way 
from a conjecture by Wigner, whose validity, in turn, is 
verified numerically for the three symmetry classes: or- 
thogonal, unitary and symplectic. The resulting wp (r) 
compares very well with the results of numerical simu- 
lations, for the three classes. The statistical analysis of 
the admittance of mesoscopic conductors provides addi- 
tional information on such systems not contained in the 
investigation of dc transport properties, and thus points 
to an interesting avenue of future research. 
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FIG. 1. Mesoscopic capacitor: A cavity (thick line) is con- 
nected via a perfect lead to reservoir 1 and capacitively cou- 
pled to a macroscopic backgate (thin line) connected to reser- 
voir 2. The cavity is ballistic and its classical dynamics is 
chaotic. 

FIG. 2. The probability density of a -the ratio of the elec- 
trochemical to the geometric capacitance- for the orthogonal 
case (Eq. (8)), for a number of values of r\. 

FIG. 3. The invariant derivative h(E). The E\ are the 
poles of K(E) and the T\ the corresponding widths. Replicas 
of K with the same width distribution (according to Wigner's 
conjecture) are generated via the transformation (12) and 
used to subdivide the area between successive levels into n 
strips. 

FIG. 4. The distribution Wp(r) of time delays for one chan- 
nel and in the absence of direct processes, for the (a) orthog- 
onal, (b) unitary and (c) symplectic universality classes. The 
dotted curves are proportional to the theoretical probability 
density given by Eq. (7). The points with the finite-sample 
error bar are the results of the numerical simulation described 
in the text: 200-dimensional matrices were used in the three 
cases. The agreement is excellent. 
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